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ABSTRACT 
> 

O 

^ ' We consider the modified restricted three body problem with power-law den- 

lO ' sity profile of disk, which rotates around the center of mass of the system with 

t~^ ■ perturbed mean motion. Using analytical and numerical methods we have found 

equilibrium points and examined their linear stability. We have also found the 

zero velocity surfaces for the present model. In addition to five equilibrium points 

there is a new equilibrium point on the line joining the two primaries. It is found 

^ ' that L2 and L3 are stable for some values of inner and outer radius of the disk 

c^ ■ while collinear points are unstable, but L4 is conditionally stable for mass ratio 

less than that of Routh's critical value. Lastly we have obtained the effects of 

radiation pressure, oblateness and mass of the disk. 

Subject headings: Photogravitational: Oblateness: RTBP:Chermnykh-Like prob- 
lem. 



1. Introduction 

The problem, after imposing a restriction as one body of the three body problem is of 
an infinitesimal(negligible) mass and remaining other two are of finite masses, is known as 
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restricted three body probleni(RTBP). The governing force of motion of the RTBP is mainly 
the gravitational forces exerted by the finite masses also known as primaries. In the RTBP 
if we take bigger primary as a source of radiation then problem called as photogravitational 
RTBP which is generalized by taking smaller primary as an oblate spheroid. 



The Chermnykh-like problem which was first time studied by IChermnykhl (119871 ). deals 
the motion of an infinitesimal mass in the orbital plane of a disk which rotates aroun d 
the center of mass of the primaries with constant angular velocity n. iGozdziewskil ( 1l998l ). 
examined the problem in the sense of nonlinear stability of equilibrium points and also 
obtained the range of parameter for the same. 

The C hermnykh-Like problem has a number of applicatio ns in different areas as celestia l 
mechanics (JGozdziewski and Maciejewskil (jl999[ ) ) , chemistry (jStrand and Reinhardtl (Il979[ ) ) 
etc. Also, the importance of the problem have seen in the extra solar planetary system 
( jRivera and Lissauerl ( I2OOOI ). Ijiang and Ipl (I2OOII ) etc.). 



Further the effect of disc is very helpful in t he st ud y of resonance cap ture of Kuiper 
Belt Objects (KBOs) as given in IJiang and Yehl (|2004J ). iPapadakisI ( l2005al ). by taking as- 
sumptions as constant mass parameter and variable arigular ve locity parameter, analyzed the 
equilibriur n point and zero velo city curve; IPapadakisI (l2005bf ) also studied the problem nu- 
merically. iJiang and Yehl (|2006[ ) examined the Chermnykh problem with /i = 0.5 and shown 
a deviation in the result of clas sical RTBP; also they have found the new equilibrium points 
in spite of Lagrangian points. lYeh and Jiang! ( 120061 ) have found the condition of existence 
of new equilibrium points analytically and numerically. 



Ishwar and Kushyahl (120061) examined the linear stability of triangular points with P-R 



drag. Again iKushvahl (l2008l ) examined the stability of coUinear points and found unstable 
points. 

Motivating by the importance and applications of the Chermnykh-like problem, we 
modeled a generalized photogravitational Chermnykh-like problem (section-|2]) in which we 
consider the angular velocity parameter n > 1 which depends on the gravitational force 
of the disk assumed, radiation force of radiating (bigger) primary and oblateness factor of 
smaller primary which is an oblate spheroid. Further we determine the equilibrium points 
(section-|3]) and zero velocity surface (section-H]) and then find the stability of the equilibrium 
points (section-|5]) finally conclude the results (section-|6]) . 
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2. Mathematical Formulation of Model 

Let us consider the motion of an infinitesimal mass governed by the gravitational force 
from the radiating body of mass mi, oblate spheroid of mass m2 [mi > 1712) and a disk 
which considered, around the central binary system, of thickness h ^ 10~^ with power-law 
density profile p(r) = ^, where p is natural number (here we take p = 3) and c is a constant 
determined by the help of total mass of the disk. 

We have chosen the unit of mass such that G(mi + 7712) = 1; the primaries are separated 
by unit distance so that unit of time obtained by the choice taken; u = — '^ — be the mass 
parameter then we have the masses Grrii = fi and Gm2 = 1 — /i. 

Let us suppose Oxyz be the rotating coordinate system having origin O at the center of 
mass of the primaries which is fixed with respect to the inertial system and angular velocity 
oj along z-axis. P(a;, y,0), 74(xi,0,0) = (— /x, 0,0) and S(x2,0,0) = (1 — /x, 0,0) be the 
positions of infinitesimal body, radiating and oblate primaries respectively, relative to the 
rotating system. 

So, the equations of motion of the i nfinitesimal mas s in x?/-plane with respect to as- 



sumptions taken above be written as[as in iKushvahl (I2OIII )] 



where 



x-2ny = fia;, (1) 

y + 2nx = Qy, (2) 



„ _ 2 (1 -/i)gi(x + /i) jj,{x + jj,-l) 3nA2{x + n~l) 



and 



r^ _ 2 _ (l-/^)gil/ _ /fy _ 3 /1^21/ _ 

\Ly — n y 3 3 5 Vy, 

'1 '2 ^ ' 2 



^^nV+l!+(W^+/i + M2_^^ (3) 

2 ri r2 2r2 



with ri = ^J{x + /i)2 + 1/2, r2 = \J{x + ^ — ly + y'^ and r = a/xM-^^. 

In above expression gi = (1 — /3) = (1 — -^), the mass reduction factor of radiating 
body, Fp is the radiation pressure force and F„ is the gravitational force of same primary. 






A2 = ^^2 ^ is the oblateness coefficient of second primary [as in lMcCuskevI (Il963f )]. Re and R 
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are the equatorial and polar radii of the same body respectively and R is the distance between 



primaries. Then the mean motion of the primaries is given as n = w gi + 1^2 — 2/(,(r) which 
is greater than 1. The potential {V) and gravitational force /^(r), of the disk are given as: 



(6 - a) 1 , 7 , log ■ 
ao r 8 r^ 



V = -2ch7i^-f^- + LchTiA^, (4) 



Mr) = -2c,J^\-\cM'^, (5) 

ab T^ 8 r^ 

where, a and b are inner and outer radii of the disk respectively. We assume that the 
gravitational force /^(r) is radially symmetric, so we have -jhif) and -fb{f) as x and y 
components of the force /b(r) respectively. Now, with the help of equations ([1]) and ([2]) we 
find the energy integral of the problem, 

C = -x'^-f + 2n, (6) 

where constant C is known as Jacobi constant. 



3. Equilibrium Points 



The coordinates of equilibrium points of the problem are obtained[as in lMoultonI (119141 )] 
by equating R.H.S.of the equations ([1]) and ([2]) both, to zero i.e. Q^ = fly = and solving 
them for x and y. In other words 

2 (1 — /i)gi(x + /i) /i(x + /i — 1) 3 fiA2(x + fi — 1) 
n X 



r\ rf 2 rf 

^ , (6 - a) X 3 , , 6 X ^ ,_, 

^2 (l-^)gii/ ^ly S fiA2y Oc/tTT ''^ ~ ^^ ^ 
rf rf 2 rf a6 r^ 

^c/^7rlog-4 = (8) 

8 ar* 

Solving equations (ITj) and ([H]) for x and ?/ we get the equilibrium points, separately given in 
the following subsections. 
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3.1. Collinear Equilibrium Points 



For the collinear points we have y = 0, so that ri = |x + /i|, r2 = |x + /i — 1| and r = |a:|. 
Now suppose 



[1 - fi)qi{x + ji) /x(x + /i-l) 3/iA2(x + /i - 1) 



f{x,y) = r?x 



, {h — a) X 3 , , h X 
—2chT[ : — -CAiTrlog — J = 0, 



ah r"^ '6 a r 

g{x,y) = n'^y 



1 - fi)qiy fiy 3 fiA2y ih - a) y 

Zcnn ; - 

ao r-^ 



rf 2 Tg 



— -c/ivrlog — 7 = 0. 



a r^ 



(9) 



(10) 



Consequently by substituting y = in equation ([9]), we have /(x, 0) = = K{x) (say) i.e. 



i^(x) = n X 



2 v^ ~ A^)Q'i(^ + /^) /i(a; + /i — 1) 3 /iA2(x + /i — 1) 



p + /i| 



(6 — a) X 3 6 X 

—zcnn : — -; — — mvr log — -; — r 



ab \x\ 



|x + /i — ip 2 |x + /i — 1|^ 
= 0. 



a \x\ 



(11) 



Now, for the sake of simplicity we first divide the plane of motion Oxy into three parts 
relative to the primaries as 1 — ;U < x , — /x < x < 1 — /i and x < — /i, we further divide 
second part into two sub parts asO<x<l — /i, — /i < x < 0, and for each considerable 
interval we have the function K{x), which we shall use in equation (TTT]) for further analysis, 
as follows: 



K{x) 



n^x 






„2 _ (l-/')gi 

ao x'^ 
(l-A')gi 



n^x 



+2c/i7r 



ih-a) 1 



+ 






M 


3 

2G 


M2 


(x+/.-l)2 


f^+Z^-l)'* 


|c/i7rlog 


6 1 






1 3 M2 
-^ 2 (x+/.-l)4 


|c/i7rlog 


6 1 






1 3 M2 

"^ 2 (x+^-l)4 


|c/i7rlog 


6 J_ 




/^ 


1 3 


MA2 



-2chTT 



(b-a) 1 



2 T^ 9 C^j_,,_1'l4 
6 1 



(a;+^-l)2 ' 2 (x+^-l) 

|c/i7rlog 



li 1 — fi < X, 
If < X < 1 - /x. 
If - /i < X < 0, 
If X < —fi, 



:i2i 



Case(l) when 1 — /x < x: Let the distance between /i and equilibrium point on the 
X-axis in the interval [1 — fi, +00) be p then x — X2 = {x + p — 1) = p > and so in this 
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case X — xi = (a; + /i) = (1 + p) > and x = (1 + p — /i) > 0. Substituting these values in 
equation (ITTl) and simplifying, we have 

8nV^° + Bip^ + B2p^ + (bs- 16chn ^ ~^' - 8p ) p^ + {B4 - ASchn ^ ~^' + 

\ ab J ab 

8 f 2c/i7r-^-^ - 5 I u + 2Ai?}p^ + {fig - ASchix ^ ~ "^ + 16 f 2c/i7r-^-^ - 5 ) u + 
\ ao / ab \ ab J 



96u^ - 2 VIp^ + {^e - 16c/i7r^ — -^ + 16 ch-n^ — -^ - 5 u + 14V - 72u^ + 8pnp 

ab \ ab J 

+B-jp^ + B^p^ + Egp + 5io = 0. (13) 

where, Bi = 48^^ - 32nV, ^2 = 120n2 - 160n2p + 48n2p^ 
53 = 160^2 - 8gi + 8 (-40^2 + gi) p + 192^^^ - 32r^V^ 
Bi = 120^2 - 3c/i7rlog ^ - 24gi + 16 (3gi - 20^^) p + 24 (-gi + 12^2) p^ - 96p3 + 8r^2p^ 



^5 = 48^2 - Gchirlog - - 24gi - 4 (40n2 + 3^2 - 18gi) p + 24 (8^^ - 3gi) p^ + 

(Ji 

24 (-4n2 + gi) p^ + 16nV^ 

6 
^6 = 8n2 - Schnhg - - 8gi - 4 (8^^ + 15^2 - 8gi) p + 12 {An^ + 3A2 - 4gi) p^ + 

32 {-n^ + gi) p=^ + 8 {n^ - gi) p^ 

Bj = -40 (1 + 3v42) p + 48 (2 + 3^2) p^ - 36 (2 + A2) p^' + 16p^ 

^8 = -8(1 + I5A2) p + 24 (1 + 9A2) p2 - 12 (2 + 9A2) p3 + 4 (2 + 3A2) p^ 

Bq = -6OA2P + 144A2p2 - 108A2p^ + 24A2p^ and B^ = -I2A2P + 36A2p^ - 36A2p^ + 
l2A2p\ 

Suppose all the quantities are constants except p on which the roots of the above 
equation depend. So, here we assume that L.H.S. of the equation flT3l) is the function of p 
and p. For p = 0, we have 

p^ {8n^p^ + ASn^p^ + UOn^p^ + C,p^ + C2P' + C3P + C4) = 0, (14) 

where, Ci = 160^2 - I6c/i7r^ _ gg^, C2 = 120^2 - 48c/i7rM _ Schnlog J - 24gi, 

C3 = 48^2 - ASchn^ - 6c/i7rlog ^ - 24gi and C^ = 8n^ - 16chn^ - 3c/i7rlog ^ - 8gi. 
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One can see that the equation flMj) has four roots equal to zero and others come from 
remaining factor. So applying the theory of solution of the algebraic equation assuming as 
/i is very small, the four roots of the equation (IH|) are expressible as a power series in jjl^^^ 
and vanish with /i. Two of the four roots are real and others are complex with the real value 
of /i^/"^. There fore the series is given as: 

12 3 4 

p = ai/i* + 02/^"' + «3Ai* + a4/W + . . . , (15) 

where ai, 02, 03, «4, • • • are constant to be determined. Putting this p into the equation 
f fT3|) and equating, the coefficient of corresponding power of p^^^, to the zero we get 



«i = ±^, «2 = 4t«S «3 = ± '^ 6^;^2 ^' «^ and 
«4 = _Mffe±g^|pM«|, where 



a 



12aM2 J _ o^A^2 



p, rfi = 8a6ra2 - 16c/i7r(6 - a) - 3abch7r\og f - 8a6gi, 



^2 = 8a6n^ + 32chn{b — a) + 9abchnlog - + IGabqi, 

ds = [AA8a%'^n^ - 6656aVc/i7r(6 - a) + 1024c^/iV^(6 - a)^ + 135a^6^c^/i^7r^(log -)^ + 
144a6c/i{15a6n^ + ichn^b - a)}log -], 

(Ji 

di = 32ab{104abn'^ + 32ch7f{b - a) + dabqichnlog f } + 256a%'^qf, 

rfs = [5120^62^^ - 16384aVc/i7r(6 - a) + 20480^^0^/1^^(6 - a)^ + 5952a%'^n^ch7r\og- + 

2232a%\'^c^h^7r^{\og -f - 27a%'^c^h^7c%log -f + 12288abn'^ c^ h\^ {b - a)log -], 

a a a 

de = 2048a6n2{4a6n2 + 10bchn{b - a) + 3abchTr\og ^}gi + 5120a%Vqf. 

So after calculating these constants and with help of equation ( TT5l) we get p. Therefore 
in the case (1) we have 

111 1 i 2 3 4 

ri = \x + p\ = l+p = l + aip4, + a;2/U4 + 03/14 + a^pi + . . . , 



^2 = \x + p — 1\ = p = aipi + a2p^ + a3/U4 + 04/14 + . . . , (16) 

111 1234 

= \x\ = \ — p + p=\— p + aipi + a2P* + cn^P'^ + a4/x4 + . . . . 



Case(2) when —p < x < 1 — p: sub case(i) when < x < 1 — p: let us suppose that 
the distance between p and equilibrium point on the x-axis in the interval [0,1 — p) be —p 



then X — X2 = {x + fi — 1) = —p < and so in this case x — xi = [x + fi) = [1 — p) > and 
X = {1— p — p) > 0. Putting these values in equation ( JTTh for same interval and simplifying, 
we have changed form of equation flT^ as follows: 



8nV'° - 5ip9 + B2p' + (-B^ + I6ch7c^^-^ - 8p) p' + {B^ - ASchn^^-^ + 

\ ah J ah 

8 (iQch-K^^^ + h\p- 2V}p6 + {-^5 + 48c/i7r^^^ - 16 f 2c/i7r^^^ + h\p + 
\ ah J ah \ ah J 

96u^ - 2V}p^ + {Bq - IQch-n ^ ~ "^ + 16 ( 16c/i7r-^-^ - 5 ) /i - 14V + 72/i^ - V}p^ 

ao \ ah J 

-Brp' + Sgp^ _ Bgp + 5io = 0, (17) 

also changed form of equation ( IT^ is 

p^ (8nV^ - 48nV^ + 120nV^ - C^p^ + C2P^ - C^p + C^) = 0, (18) 

where Bi,B2, . . ■ ,Bg, Bio, C*i, C*2, Cs and C4 are given as in case (1). Similar analysis as in 
case (1) provides 

111 1 i 2 3 4 

ri = |a:; + p| = l — p = l — Q;ip4 — 0:2^* — asP* — a4p4 — . . . , 

I 1 I 12 3 4 /I n\ 

r2 = \x + fi — 1\ = —p = —aipi — a2p,* — a^p,^, — a^p* ~ • • • , U^J 

111 1 1234 

r = |x| = l — p — p=l — p — aipi — a2P^ — ctsP* — a^pi — . . . . 
where constant ai, 0^2, as, 04, . . . are given as: 

«i - ±^^, a2 - 4^(-ajS as - ± ^x^^ (-«)* and 

04 = — ^^^ — ^ ^ 12A d^ ^ — —{—a)^ and a, di, ^2, ^3, ^4, (is and d^ are as in case (1). 

Sub case(ii) when — p < x < 0: let the distance between p and equilibrium point on the 
X-axis in the interval (— p, 0) be — p then x — X2 = (x + p — 1) = — p < 0, x — xi = (x + p) = 
(1 — p) > and x= (1 — p — p) <0 (The main difference in sub case (i) and(ii) is that 
in (i) X > and in (ii) x < 0). Putting these values in equation flTTj) for same interval and 
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simplifying, we have changed form of equation flT3|l of the case (1) as follows: 

8nV'° - B,p^ + 52p« + (-B-, - Wchn^-^^ - 8fi] p' + {B, + ASchJ^^ - 

\ ah J ah 

8 (2c/j-7r '^ ~/^ - 5 1 u - 2V}p^ + {-B^ - 48c/i7r-^-^ + 16 f 2c/i7r-^^^ - 5 ) u + 



96u^ - 24unp^ + {^e + 16c/i7r^ — -^ - 16 ch-n^ — -^ - 5 u - 14V + 72/i'^ - 8/i^lp 

ao \ ah J 

Brp' + Bsp^-Bgp + B,o = 0, (20) 

also, changed form of equation (fT4|) of the case (1) is 



/ 



8nV' - 48nV' + 120^^' + (Ci + 16gi) p' + ( C2 + 96c/i7r^^-^ ) p^ 



C3 - Uchirlog - + 48gi ) p + ( C4 + 32c/i7r *^ "'* 



, .„,, , , , , ^. , , , 0, (21) 

a J \ ah 

where Bi, B2, . . . , Bg, Biq, Ci, C2, C3 and C4 are given as in case (1). 

- - = - 

ri = |x + p| = l — p=l — aip4 — a2p4 — a3p4 — a4p4 — . . . , 

r2 = |x + p — 1| = — p = — aip4 — a;2p4 — ttsP* — a4p4 — . . . , (22) 

/I \ 1 12 3 4 

r = \x\ = —{I — p — pj = — 1 + P + aip4 + Q;2p4 + a3p4 + a4p4 + . . . . 

where constant ai, 02, as, 04, . . . and a given as in case (1) whereas 
di = Sahn^ + 16chn{h — a) — ^ahchnlog - — Sahqi, 
d2 = Sahn^ — 32ach'n-(h — a) + dahchiilog - + 16a6gi, 

h 
4 = [448a^6^n^ - QQ^Qahchn^Trih - a) + 1024c^/iV2(6 - a)^ + 135a%'^c'^h'^'K^{\og -^ + 

h 
lA'iabchllbab'n? - Ach'K{h - a)}\og -], 

Ci 



4 = 32a6{104a6n2 - 32c/i7r(6 - a) + 9ahchiT\og ^}qi + 25Qa%'^ql 



4 = [512a%\^ - 16384a6n^c/i7r(6 - a) + 20A80n^ c^h\^{h - af + 5952a%'^n^ch7rlog - + 

(Jj 

2232a%^n^c^h^7i\\og -f + 27a%\^h'Ti\\og -f - 12288ahn^c'hV{h - a)\og -], 

a a a 

21,2^2^2 



de = 2048ahn^{Aahn^ - 10c/i7r(6 - a) + 3a6c/i7rlog ^}qi + 5l20a%^n^ql 
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Case(3) when x < —fi: let the distance between 1 — /i and equihbrium point on the 
X-axis in the interval (—00, — /i) be denoted by 1 — p then x — Xi = {x + ft) = {p — 1) < is 
negative and so x — X2 = (x + /i — 1) = (p — 2) < and x = (p — p — 1) < 0. Putting these 
values in equation fITT]) and simplifying we get 

8n2p'° + Bip^ + B^p^ + 53p^ + 54p' + 55p^ + B^p^ + Bjp^ + B^p" + B^p + 5io =(0?) 

where in this case, 

Bi = -IVlr? - 32n^fi, B2 = 696^2 + 416n^fi + 48n2p^ 

53 = -2528^2 + I6c/i7r^ + 8gi - 8 (296^2 + gi - 1), 

54 = 5944^2 - 176ch7c ^ ~^' - Schirlog - - 88gi + 8{968n2 - IGch-K ^ ~ "^ 

ab a ab 

-8gi - 9}p + 24 (I24n2 + gi - l) p^ + 352nV^ + 8n2p^ 

55 = -9456^2 + SlQchir ^ ~^' + SOchirlog - + 408gi + 41-4008^2 + AOchn ^ ~^' + 

ab a ab 

3A2 - 42gi + 68}p - 8 {lOSOn'^ + 27qi - 24) p^ - 24 {68n^ + gi - l) p^ - 80n2p^ 

Bq = 10312^2 - 2064c/i7r ^ ~/^ - 123c/i7rlog- - 1032gi + 4(5448^2 - IQAchn ^ ~/^ 

ao a ao 

-15^2 + 12gi - 140}p + 12 (1284^2 - 3^2 - 64gi - 52) p^ + 8 (516^^ + 26gi - 2l) p^ 
+8 (41^2 + gi - 1) p^ 

B7 = -761Qn^ + 3072ch7i ^ ~/^ + 264c/i7rlog - + 1536gi + 8{2432n2 + UGchn ^ ~/^ 

ab a ab 

+I5A2 + 72gi + 85}p + 4 (-4320^2 + 36^2 + 336gi + 264) p^ + 4 (-1536^2 + 9A2- 
176gi + 114) p^ - 16 (44n2 + 4gi - 3) p^ 

^8 = 3648n2 - 2Q88ch7T ^ ~/^ - 312c/i7rlog - - 1344gi + 8(1376^2 - 208ch7r^-^-^ 

ab a ab 

-I5A2 - 144gi - 61}p + 8 (I488n2 - 27^2 + 144gi - 123) p^ + 4 (1344^2 - 27^2 + 288gi 
-150) p^ + 4 (208n2 - 3^2 + 48gi - 13) p^ 

Bg = -1024n2 + 1280ch7T ^ ~/^ + 192c/i7rlog - + 640gi + 4{-896n2 + 256c/i7r ^ ~^' 

ab a ab 

+ I5A2 + 224gi + 48}p + 8 (-576^^ + I8A2 - 48gi + 6O) p^ + 4 (-640^2 + 27^2 - 224gi 
+96) p^ + 8 (-64^2 + 3^2 - 32gi + 12) p^ 
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Bio = 128^2 - 256c/i7r-^ — -^ - ASchnhg - - 128gi + 4(128^2 - GAchir^ — -^ 

ab a ah 

-3^2 - 64gi - 8}/i + 12 (64^2 - 3^2 - 8) /i^ + 4 (128^^ - 9^2 + 64gi - 24) /i^ 
+4 (32n2 - 3^2 + 32gi - 8) /. 



Putting /i = in equation (!23|) . we have 

(-2 + p)^(-l + p)2[8a6nV^ - 32a6nV + 48a6nV{-32a6n2 + 16c/i7r(6 - a) 

+8a6gi}p + 8a6r2^ — 16c/i7r(6 — a) — 3a6c/i7rlog 8a6gi] = 0. (24) 

a 



Clearly out of ten roots of equation fl24j) . four equal to 2, two equal to 1 and others come 
from remaining factor. Now here we want to apply the theory of solution of the algebraic 
equation similar as earlier cases assuming as p very small and L.H.S. of equation fl2^ is 
function of p and p only. So the four roots of the equation (123|1 are expressible as a power 
series in p^^^ with an extra term as 2 and two roots as a power series in /i^/^ with an extra 
term as 1 and these series not vanish with p due to extra term. The power series in p^/"^ 
with an extra term as 1 is independent to p as its constant coefficient comes out zero, so this 
series is unimportant whenever the power series in p^/'^ with an extra term as 2 depends on 
p. The two of the four roots are real and other are complex with the real value of /x^^^.The 
power series is given as: 



12 3 4 



p = 2 + aip4 + a2^4 + 0:3^4 + 0:4^4 + . . . , (25) 

where ai, 0^2, as, 04, . . . are given as in case (1) which is obtained by putting p from fl25|) 
into the equation ( !23l) and equating, the coefficient of corresponding power of p^^^., to the 
zero. In this case a is same as in earlier cases and ^3, d^ are as in sub case (ii), whereas 

di = Sahn^ + IQchnih — a) — Sabchrclog - + Sabqi, 

(^2 = —Sabn^ + 32c/i7r(6 — a) — 9abch7rlog - + IGabqi, 

d4 = -32ab{104:abn^ - 32c/i7r(6 - a) + dabchnlog ^}qi + 256a%'^ql 

de = -2048abn^{4:abn^ - 10chn{b - a) + Sabchnlog f }gi + 5120a%'^n^qf. 

Again from equation fl2Sl) we get p and there fore in this case we have 

2 3 4 

a2p^ — ot-iP'^ — a^pi — . . . , 
aipi — a2P^ — Q-sP^ — a^pi — ■ ■ ■ , (26) 

12 3 4 

4 — Q;2p4 — 0:3^4 — 0:4^4 — . . . . 



n = 


■■ \x + p = - 


-(p- 


-l)- 


= -l- 


-aip4 


r2 = 


-■ \x + p— l\ 


= - 


-{p- 


-2)=- 


-aip-i 


r = 


\x = -(p- 


- p- 


-1) = 


=-l+ 


p — a 
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For numerical calculation we set the parametric values as follows: Mass parameter(/i) 
_ = 0.000953728, where M,(Sun mass)= 1.98892 x 10^°^^ and mj(Jupiter mass) 



l.mS7xW^kg, radiation factor (gi) = 0.75, oblateness effect(A2) = 0.0025, disk's mass(M6) = 
0.4 (taking disk's inner radius(a) = 1, disk's outer radius(6) = 1.5, control factor of density 
profile(c) = 1910.83, disk's thickness(/i) = 0.0001 and vr = 3.14). Putting these values into 
the equation (ITTl) for each intervals and solving we have collinear equilibrium points with 
respect to the relevant intervals given below: 

In interval (1 — /i, +cxd) we have six real values of x out of that five lie outside the interval 
and remaining one i.e. Li = 1.05667 belongs to the interval. Similarly in case of interval 
(0, 1 — yu) we have four real values of x, one of them i.e.L2 = 0.813609 belongs to the interval 
(0, 1 — /u) and three other lie outside the interval. In the interval, (— //, 0) we again have 
four real values of x, one of them i.e./i = — 0.00879106(new point) belongs to the interval 
(— /i,0) and three other lie outside the interval but in the interval (— oo, — /i), we have only 
two real value of x one of them i.e. L3 = —0.823420 belongs to the (—00, — /i) and other one 
is positive. 

We also find the collinear equilibrium points with the help of equations f lT6|) . flT9l) . fl22|) 
and fl26|) . after evaluating ai,i = 1, 2, 3, 4 . . . and hence p for each case, also for mean motion 
(n) we set here disk's reference radius r = 0.99, given as Li = 1.04411, L2 = 0.971566, 
/i = —0.975463 and L3 = —0.977926 which is very close to the points obtained above one 
excepting the new point /i which is also outside of the relevant interval. 

In figure [1] (a), /(x,0) = K{x) Vs x contains four curves (i), (ii), (iii) and (iv) for the 
different intervals of x or domains of /(x, 0) but same range of /(x, 0). These four curves 
intersect x-axis at one and only one point i.e. at Li, L2, L^ and /i in their respective intervals 
(but out side the respective intervals we find some more intersections points of each curves 





Fig. 1. — Collinear Points: /(x,0) Vs x- at gi = 0.75, A2 = 0.0025, a = 1, 6 = 1.5 (a) range 
and domain of /(x, 0) both different (b)range and domain of /(x, 0) both same. 
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X 



Fig. 2. — Comparison of nature of curve /(x, 0) drown at different values of parameters: 
gi, A2 and b. 
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which are other roots of the flT3l) . fITTj) . fl20|) and fl23|l . These intersection points i.e. Li, L2, 
L3 and li are much clear in figure [1] (6) which is drown at same domain and range, which is 
nothing but the respective intervals of each curves. In figure [H we plots same curves as in 
figured] (a) but at different values of disk's outer radius b i.e at 6 = 1.0, 1.2, 1.5, and 2.0 in 
addition with classical case (gi = 1, ^2 = and a = b = 1.0). In other words figure [21 shows 
the nature of curve f{x, 0), at different widths of disk, which is similar in nature as in figure 
[1] (a) and intersects a;- axis at only one point Li, L2, L^ and /i in their respective intervals 
while in classical case (even in our cases also i.e. if gi = 0.75, A2 = 0.0025 and a = b = 1) 
the curves f{x,0) intersect x-axis at three different points Li, L2, L3 and for x in (— /i,0), 
/(x, 0) does not intersect a;- axis. 

Thus, we conclude that in the presence of disk, there are always four different intersection 
points Li, L2, L3 and /i of curve f{x, 0) with x-axis, each lies in the four different intervals 
described above. 



3.2. Triangular Equilibrium Point 



In this case we have y ^ 0, for the convenience let us suppose that ri = g^ (1 



and r2 



1 + 62, where 61,62 <^ 1. Putting these values of ri and r2 into the equations 



X + yu)^ + y"^ and r| = (x + /i — 1)^ + y"^ and then solving with the rejection of second 



and higher order terms of ^i and 62, we get 



2/3 



X 



2/3 



^2 /" + (^1 ^1 ^ '^2) and y = ±q^ 



1/3 



2/3 



Qi 



2/3\ 



)6i + 62 1/2. We de- 



termine the values of 61 and 62 by putting the values of x, y (after neglecting the terms 
containing 61,62 <^ 1), ri = q/ (1 + 61) and r2 = 1 + ^2 into the equation ([7]) and ([8]) and 
solving them with the rejection of second and higher order terms of 61 and 62, we get 

(b-a) 1 , 3 , log' 



61 



n 



2ch 



11- 



ab 



-chn- 



ifi' + qr {I -fi))3/2 8 (/.2 + gP(l_^))2 



2/3/ 



,(27) 



3(1 + 1^2) 



3 . 2 ^ , (b — a) 

2 ab 



1 



(/i^ + gf(l-/i))3/2 



+ 



-ch 



11- 



log- 



8 (^2 + ^2/3^^_^^^2_ 

Hence, the coordinates of triangular points are 



(28) 



X 



2/3 



- A* + (^l^^*^! - ^2 



y = ±^1 



1/3 



1 - C + (2 - qT)6, + 62 ] '/' 
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where Si and S2 are given by equation fl271) and (l28l) respectively. 

Numerically the co-ordinates of triangular equilibrium points are L4 = (0.366171, 0.641213), 
and L5 = (0.366171, —0.641213) which are obtained by using same parametric values as in 
collinear case into the equation (Q and ( ITOj) and solving them for x and y. We also calculated 
L4 = (0.447217,0.515281) and L4 = (0.447217, -0.515281) with the help of equation ([29]) by 
evaluating 61 and 62 from the equations ( 127|) and (128|) . The mean motion (ri) is calculated 
at disk's reference radius (r) = 0.99 and found very similar results. We have seen that the 
1/4(5) cire no longer remain triangular equilibrium points as they are in classical case. 



4. Zero Velocity Surface 

Since, equation ([6]) generally denotes the relation between coordinate and velocity of 
infinitesimal mass with respect to the rotating system, so by taking velocity term as zero in 
equation (E]) we have the equation of zero velocity surfaces as follows: 



2n = c 



(29) 



that is 



x^ + y^- "^i^^ 



ri 



2£t 

r-2 



HA2 



+ Ach'K 



(b-a) 1 _ 7 
ab r 4 



ri = a/(x + jj,y + y'^ + z^ 
r2 = ^/{x + jj, - 1)2 + 1/2 + 
r 



A/x^ + y^TTz^ 



chn 



log 



C 



Here, we try to know the approximate form of zero velocity surface [as in iMoultonI ( 119141 )] by 
analyzing the shape of the curves obtained by intersection of surface fl30|) with the xy-plane. 
The equation of curve is given by putting 2; = in fl30|) which is as follows: 



2 2 2(1 — a)qi 

x^ + y^+ ,r = + 



2// 



+ 



flA2 



Achn 



v/(x + /i)2+|/2 ^(x + /i - 1)2 + y2 ((a; + ^_ 1)2 + ^2)3/2 

a6 ^a;2 + y2 4 ^2 + 1/^ 



+ 



C 



(30) 



Here if we increases the value of x and y, all terms, except first and second, on L.H.S. of 
(l30|l . become negligible compare to first and second terms. Thus for the large value of x and 
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c 




Fig. 3. — Zero velocity surface: gi = 0.75, A2 = 0.0025, a = 1, 6 = 1.5 and r = 0.99. 
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y, curves become very close to a circle of radius a/C — e having equation 

2(1 - /i)gi 2/i /iAs 



x^ + y^^c 



v/(x + ^)2 + t/2 ^(x + ^-l)2+y2 ((a; + ^- 1)2+2/2)3/2 



Achn 



(b-a) 



+ Icte-'°«^ 



a& v/x2 + 1/2 4 x2 + y 



C 



(31) 



where, e is very very small. For the small value of x and y, the terms x'^ + y'^ will be negligible 
compare to remaining terms on L.H.S. of (1301) and therefore, the curves become very close 
to an equipotential surface given bellow 



2(1 -^)gi 



Achir 



{b — a 



+ 
1 



2/i 



I^A 



V/(X + /i)2 + y2 ^(x + /i - 1)2 + y2 ((a; + /i- 1)2 +^2)3/2 

b 



—Qfl-Jl 2 

ab ^x2 + 1/2 4 x2 + 2/2 



C-x2-t/2 



c 



(32) 



where e is very very small. 

The surface is given in figure |3] where two singular regions are shown by conic shapes. 



5. Stability of Equilibrium Points 

For the stability of equilibrium points (xe, ye) let us assume a small change in its coor- 
dinate as X = Xe + i,y = ye + Vy where the displacements ^ = Pie'*'*, r/ = P2e^* are very small. 
Pi, P2 and A are parameters to be determined. Putting these coordin ates into equations (HI) 
and (p) , we have the transformed equations of motion as follows [as in [Murray and Dermott 
J2000f )] 



^-2nrj = ^nl + V^l 



xyi 



f) + 2< = ^nl + r^n^, 



(33) 
(34) 



where superfix denotes the corresponding value at equilibrium point. Now putting ^ 
Pie^^^T] = P2e^^ and simplifying we have 



(A^ - fiL)Pi + {-2n\ - fi°^)P2 = 0, 
{2n\ - fi;,)Pi + (A^ - VPJP^ = 0, 



(35) 
(36) 



for nontrivial solution we have 



2n\ - nl 



-2n\ - n% 



0. 
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Simplifying above we have characteristic equation as follows: 



Now 



r^ Tg 2 r2 ao r'^ 

3 , , 6 1 3gi(l -u)(x + u)2 3u(x + u-l)2 15 uAsfx + u - 1)^ 

, {h-a)x^ 3 , , 6x^ ,„„, 

^ _^2 gi(l-^) /^ 3/iv42 2c/i7r '^^ ~ ^^ ^ ^c/irrlo -— + 

yy ^3 j,3 2 rf ah r^ 8 ar^ 

3gi(l-^ 3/ii/2 15/xA2y2 (&-a) y^ , 3 6 2^ 



, ^ r| 2 r2 a6 r^ ' 2 ^ a r^ 

_ _ 3gi(l -/i)(x + /i)y 3/i(x + /i-l)y 15 /iA2(x + /i- l)y 



r'l Tg / r2 



^xx = n -\ 5 1 — T H J h 4c/i7r 



O /ytO /y*^J 



ar 



(39) 



(6-a)x|/ 3 6x1/ 
+6c/i7r ; + -c/ivrlog -. (40) 

In the following subsection we examine the stability of equilibrium points. 



5.1. Stability of Collinear Equilibrium Points 

For the stability of collinear points {y = 0), we have 

2 , 2gi(l-;u) ^ 2/i ^ 6/iA2 , ,^,i^^ib-a) 1 



rj" r2 Tg ab r^ 

9 6 1 

+ c/lTTlog--, (41) 

8 ar^ 

_ 2 gl(l-/i) /^ 3/1^2 0,^ (^-^) 1 

r^ Tj 2 Tg ao r'^ 

3 , , 6 1 
--chTilog--, (42) 



also we have ri = |a; + /i|, r2 = |x + yU — 1| and r = |x|. To insure the stability of equilibrium 
points, we must have ^ = Pie^^,-!] = P2e^^ which can be expressed in a periodic functions. 
In other words the four roots Ai, A2, A3 and A4 of the characteristic equation fl37j) must be 
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pure iinaginary[as in iBoccaletti and Pucaccd (Il996[ )] otherwise we have an unstable point. 
Hence for this we examine here nature of roots which is easily done by finding the sign of 



r2° and fi,,,, as follows: 



D 



In case(l) when 1 — /i < x, we have fl^^ > and Qyy < and hence the discriminant 
(4n^ — i7°^ — ^yy)'^ — 4:Q'l^Qyy has positive sign. Therefore in this case, the two of the 

four roots Ai, A2, A3 and A4 of the characteristic equation ( 137|) are complex conjugate pair 

and remaining two are pure imaginary conjugate pair. That is 



Ai — —As 



\ 



f4n2 - f]0 - fiM + a/ f4n2 - fio - n^ ) - 4^° f]o 

V XX yy J V V XX yy / xx y 



yy 



is a complex conjugate pair and 



As — — A4 



\ 



(4n2 - n», - ni,) - J(4n= - a»„ - il«J' - 4!);.fiJ, 



is a pure imaginary conjugate pair. Thus the condition of stability failed. Similarly, analyze 
the case (2)(i) when 0<x<l — /iwe have fi^^ > and Qyy > 0, in case (2)(ii) when, 
— /i < x < we have Q^^ > and flyy < and in case (3) when, x < —fi we have fi'^^ > 
and flyy > 0. That is in each case we have a positive discriminant. Therefore, we conclude 
all the coUinear equilibrium points are unstable. We examined the stability of coUinear 
equilibrium points with respect to same parametric values setting above and conclude that 
generally Li, L2, L3 and /i(new point) are unstable points but L2 is stable for the disk's 
width (a = 1, 1.446 < 6 < 1.767) also L3 is stable for the disk's width (a = 1, 1 < 6 < 1.08). 



5.2. Stability of Equilibrium Point L4 

In case of L4, we have 



4^2 _ ^0 _ qO^ ^ ^2 



yy 



3fiA2 3 ^ , b 1 

— ch7clog--j, 



r. 



2,0 



^^xx^^yy ^^ xy 



9/i(l - /i)7o. 



ar, 







(44) 
(45) 



where 



7o = 2/o 



ch 



Qi 



'1,0' 2,0 






0' 2,0 



1 + 



5Ao 



2r2 
^'2,0 



n{i-fi) r 5A2 



2r2 
^'2,0 



qifi /^ fe-g log &- logo 



/yiO /yiO 



ah 



2ro 



b — a log b — log a 



ab 



2ro 



(46) 
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The superfix indicate that values are at equilibrium points. Now, from the characteristic 
equation ( [37j) we have 



(4n2 - fio^. - fioj ± ^(4^2 _ ^0^ _ QO^Y _ 4 (^o^f^o^ _ ^o^j 



(47) 



It is seen earlier that ^ = Pie , rj = P2e will be periodic and bounded if A 
pure imaginary, i.e. A^ < 0. Therefore for stable solution we must have (4n^ — f2°^ 



4^0 fio 

' XX yy 



fi°' 



In other words 



n 



chn log 



b 1 



> 36/i(l — /i)7o. For the 



1,2, 3, 4 is 

classical values i.e.(gi = 1, A2 = 0, b = a, n = 1, x = -^ — fi, y = ±^, ri = r2 = 1), the 
above inequahty reduces as 27yu(l — /i) < 1 which provides fi < 0.0385209 = /ic the value of 
critical mass in classical case. But in our case the value of critical mass for different values 
of gi, A2, and b will be obtained by following graphs. 

In figure m we plot fi Vs b, where (a) for qi = 1 and (b) for qi = 0.75 both containing 
three curves (i) A2 = 0, (ii) A2 = 0.0025, (iii) A2 = 0.0050 have same nature i.e. as we 
increases the value of disk's width b, critical mass fic also increases. Initially these curves 
increases slowly but when 6 > 1.6 they increase strictly as depicted in figur^ These shows 
that when we increase the width of disk, stability region spans slowly then rapidly after 
b = 1.6 upto fie = 0.5. 



6. Conclusion 

We have studied the dynamical properties of modified restricted three body problem 
and found that there exists a new equilibrium point in addition to five equilibrium points 
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Fig. 4.— Variation of critical mass: (i) A2 = 0, (ii) A2 = 0.0025, (iii) A2 = 0.0050. (a) at 
g^ = 1, a = 1, r = 0.99, (b) at qi = 0.75, a = 1, r = 0.99. 
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in classical problem. The zero velocity surface have obtained and examined the stability 
of equilibrium points. It is found that the L2 and L3 are stable for certain values of inner 
and outer radius of the disk and other coUinear points are unstable while L4 is conditionally 
stable upto the Routh's value of mass ratio. Thus stabihty region spans with the width of 
the disk. Hence, we conclude that the shape and size of the disk are very significant for the 
motion of the bodies in space. 
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